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QQ ' The evidence is reviewed for the mechanism of colour confinement in QCD by dual 

Q\ , superconductivity of the ground state, i.e. by condensation of monopoles. 

o\ 

_ ' 81. Introduction 

X) ■ 

<D , 

Ph . Most of the existing tests of Quantum Chromodynamics (QCD) as theory of 

strong interactions come from short distance phenomena (Deep inelastic scattering, 
e + e~ — * hadrons, Jet physics,. . . ). Perturbation theory is supposed to work in that 
regime because of asymptotic freedom. 
' On the other hand it is known that the renormalized perturbative expansion 

does not converge, not even as an asymptotic series, due to bad infrared behaviour ^ . 
Fock vacuum is unstable: quarks and gluons, which are the elementary excitations of 
£S) ' perturbation theory, never appear as asymptotic states. This phenomenon is known 

as colour confinement. 

A convincing experimental evidence for colour confinement is the upper limit on 
the cosmic abundance of relic quarks, n q : 
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n„ 



< 1(T 27 (1-1) 



n p is the abundance of nucleons. 



Eq.(lT) correspond to Millikan like analysis of ~ 10 2 gr of matter. For non 
^ . confined quarks the standard cosmological model predicts 

< 1(T 12 (1-2) 

n p 

A non perturbative formulation of the theory is needed, as well as a theoretical 
understanding of why perturbation theory works at all at short distances. 
Lattice formulation provides that formulation. 

The Feynman path integral defining the theory is regularized by discretizing 
space time, and computed numerically by Montecarlo techniques. Of course numer- 
ical computations have not the logical trasparency of mathematical derivations, but 
they can help understanding anyhow. Lattice can be used of course to compute 
observable quantities (masses, matrix elements. . . ) from first principles (Lattice for 
phenomenology), but also as a tool to explore mechanisms and structures of the 
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theory by means of "Gedanken experiments" on the configurations produced by nu- 
merical simulations (Lattice for theory). 

We will be concerned with an investigation of the second type (Lattice for the- 
ory) , and test on the lattice the idea 2 ) " 4 ) that confinement of colour is produced by 
Dual superconductivity of type II of the QCD vacuum. 

The basic idea is that the chromoelectric field produced by a QQ pair is chan- 
neled by dual Meissner effect into Abrikosov flux tubes in the same way as magnetic 
field is confined in usual superconductors of type II. 

The energy is proportional to distance 

E = aR (1-3) 

and this means confinement, a is the string tension. 

The world "dual" here means that the role of electric and magnetic quantities 
is interchanged with respect to ordinary superconductors. 

The idea is theoretically appealing in many respects. 

1) Superconductivity is a symmetry deconfining transition is a change of sym- 
metry. An order (disorder) parameter can be defined and used to explore super- 
conductivity. Dual superconductivity means a condensation of monopoles in the 
ground state. Vacuum has no definite magnetic charge, but is a superposition 
of states with different values of it 6 )' 7 ). 

A disorder parameter for dual superconductivity will then be the vacuum ex- 
pectation value (vev) of any operator carrying magnetic charge. 

2) The very existence of monopoles implies that the theory is compact, because of 
the Dirac quantization condition. Topology plays then an important role, since 
Dirac strings make the connection of the gauge field non trivial. A formulation 
in terms of parallel transport, like Wilson's lattice theory, is then superior with 
respect to the perturbative formulation in terms of local fields, which misses 
topology. 

3) Monopole charges are always £7 (1) 8 -*' 9 -*: in SU(N) gauge theory there are N —1 
£7(1) magnetic charges. Monopole species are associated to any field # in the 
adjoint representation 

$( x ) = Y J $ a \ a (1-4) 

a 

\ a are the generators of the gauge group in the fundamental representation. 
In field configurations monopoles are located at the sites where two eigenvalues 
of the N x N matrix <P, eq.(l-4), coincide. This is of course a gauge invariant 
statement. The N — I £7(1) fields associated with magnetic charges are exposed 
by a gauge transformation £7 (x) which diagonalizes <P(x). Such a transformation 
is called an abelian projection. 

There are a functional infinity of choices for and to each of them monopole 
species are associated. What species do condense in the vacuum to produce 
superconductivity is a dynamical question, and is actually what has to be in- 
vestigated to understand the mechanism of confinement. A possibility is that all 
species are equally important for confinement and condense in the vacuum 10 ) . 
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One tool to investigate this issue is to define a disorder parameter for different species 
and to measure it in connection with the deconfining transition. As we shall see below 
this connection can be made unambigous by a measurement of the critical indices 
or effective critical indices of the phase transition, extracted from the behaviour 
of the disorder parameter. A somewhat different attitude is to look at the abelian 
dominance n ) . For a given abelian projection physical quantities as the string tension 
are measured in the full theory and in the U (1) theory resulting from the abelian 
projection. If the latter determination is a good approximation to the exact quantity, 
people say that there is abelian dominance. This happens to 80% approximation in 
the so called maximal abelian projection. In addition one can separate the abelian 
field into a part due to monopoles plus a residual part with no topology 12 ) . If the 
monopole part dominates people say that there is monopole dominance. Again this 
happens in the maximal abelian projection. Monopole dominance is then considered 
as a strong indication that confinement is due to monopoles. 

The two approaches are in our opinion both important to disantangle the struc- 
ture of the theory. 

Our strategy is the following 

a) We define a disorder parameter for U(l) dual superconductivity, and test its 
construction in compact U(l) gauge theory. The costruction is also tested with 
other well known systems, like the XY 3d model, which describes the transition 
to superfluid He^. 

In both cases we are able to detect the change of symmetry of the ground state, 
and to determine the critical indices. 

b) We then define a similar disorder parameter for abelian projected monopoles of 
non abelian gauge theories, and explore by it the occurrence of dual supercon- 
ductivity in connection with confinement. 

A systematic analysis is in progress. Our preliminary results confirm that monopoles 
defined by different abelian projection condense in the confined phase, and the su- 
perconductivity disappears in the quark gluon phase, supporting the view of ref. 10 ) . 

We will start these lectures by a brief introduction to basic superconductivity in 
sect. 2. 

In sect. 3 we shall recall the main properties of monopoles and the concept of 
duality. 

We will then describe the disorder parameter for condensation of monopoles in 
U{\) gauge theory and for condensation of vortices in 3d XY model (sect.'s 4,5). 

The abelian projection and the physical meaning of monopoles will be discussed 
in sect. 6. The results for SU(2), SU (3) and the other evidences from lattice for dual 
superconductivity mechanism will be reviewed. 

Sect. 7 will summarize the state of the art and present the open problems. 
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§2. Superconductivity as a symmetry. 

A relativistic version of the Ginzburg-Landau free energy, which is the statistical 
analog of effective action, is 

C = -\f, u F, u + (D^)\D^) - V(0) (2-1) 

where ^ is a complex (charged) scalar field describing Cooper pairs of charge q = 2e. 

= (d„ - iqAJ* (2-2) 

is the covariant derivative, 

V(0) = ^0-f, 2 ) (2-3) 

the potential, with p, and A functions of the temperature T, p = p(T), A = A(T). 

Minimizing C defines the ground state. If p 2 > the minimum corresponds to 
some {<&) 0, or to the Higgs phase. At T where /j = 0a transition to normal phase 
{p 2 < 0, (#) = 0) takes place. (#) is the order parameter of superconductivity. 

Putting <P = pe te , p > 0, under gauge transformations of angle a(x) 

p(x) -> p(x) 9(x) -> 6(x) + a(x) A^(x) -»■ ^(s) - <9 M a 

Therefore for the covariant derivative 

= e ie % - iq(A^ - 8^0)] p 

A^ = (A^ — d^Q) is gauge invariant and F^ u = d^A v — d v A^, since d^d does not 
contribute. 

C can be rewritten as 

1 771 2 ~ ~ - 

C = --Fp,Fp, + —A^An + C [p] (2-4) 
and the equations of motion for the electromagnetic field read 

d\xFjiv + m 2 A v = (2-5) 

with m 2 = 2q 2 (<P) 2 . 

In a stationary state with no charges Aq = 0, doA = and equation (2-5) gives 
(H = V A A) 

V AH + m 2 A = (2-6) 
V 2 H + m 2 H = (2-7) 

Eq.(2-6) means that a permanent current (London current) j = m 2 A exists, and 
since E = and E = pj, p = 0. If m 2 , or (<2>) is different from zero, there is 
superconductivity (p = 0) 

Eq.(2-7) is nothing but Meissner effect: the penetration depth of the field H is 
A = 1/m and is again finite if (#) / 0. 
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A side consequence of Meissner effect is flux quantization: outside a flux tube, 
at distances larger than A, A = 0. The integral around a circle C centered on the 
section of the tube of A is zero 

= <j>Adx = <j>{A-V0)dx 
or, since the flux F(H) = § Adx 

F(H) = H= (2-8) 

The key parameter is the order parameter {$) which signals the Higgs phenomenon. 

(<£) / means condensation of charge. Indeed, if the ground state has a definite 
charge, the expectation value on it of any charged operator C is zero: (0|C|0) = 0, 
since C|0) belongs to a different eigenvalue of the charge than |0). 

Superconducting vacuum is indeed known to be a coherent superposition of 
states with different numbers of Cooper pairs 6 )' 7 ). 

There are two characteristic lenghts in the system: the correlation length of 
the field, or the inverse Higgs mass A = 1 /M, and the penetration depth of the 
photon, A. If A > y/2A the superconductor is called type II, and the formation of 
Abrikosov flux tubes is favoured in the process of penetrating the material with a 
magnetic field 13 ). If the opposite inequality holds, A < V2A, when the magnetic field 
is increased there is an abrupt penetration of it at some value and superconductivity 
is destroyed. 

In principle many independent charged fields could condense in the vacuum. In 
that case 

^ = E«?<*<> 2 (2-9) 
§3. Monopoles and their topology. 

3.1. U(l) monopoles. 

Maxwell's equations in the presence of both electric and magnetic currents, j^, 
3 Mi are 

d^=f d^=f M (3-1) 
F^ v is the familiar field strength tensor, F* v its dual 

If both j v and j M are zero the transformation 

Ffa, cos 9Fp, + sin 6F* U (3-2) 
F; u - cos 6F; u - sin 9F^ (3-3) 

is a symmetry of the system for any value of 6. In particular for 6 = tt the transfor- 
mation becomes 
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which is known as duality transformation. 

In nature, within the present experimental limits, j^j = 0, since no isolated 
magnetic charge has been found. Therefore 

d^F*^ = (3-4) 

the general solution of this equation is 

iV = - ( 3 - 5 ) 

which makes eq.(3-4) identically satisfied (Bianchi identity). 

Of course if j^- 7^ Bianchi identity is violated, and the vector potential 
cannot be defined. The way out of this difficulty, was proposed by Dirac 14 ). A 
monopole can be seen as the end point of a magnetic flux tube, an infinitely thin 
solenoid (Dirac string), thus preserving Bianchi identities. The solenoid will be 
physically invisible if for any particle with charge e the parallel transport around it 
is trivial 

exp(ie (p Adx) = 1 or <P ■ e = 2irn (3-6) 



The magnetic flux on the other hand is related to the visible magnetic charge M of 
the monopole as <P = M, and this requires 

eM = 2vrn (3-7) 

with n integer, for any particle. Any charge is then multiple of the same elementary 
charge 

-s < 3 ' 8 > 

Theory is compact. A formulation in terms of parallel transport is then superior in 
that it naturally describes the non trivial connection of space time produced by the 
presence of Dirac strings. 

One could formulate the theory in a dual form, by introducing a dual vector 
potential, and then dpF^ = play the role of Bianchi identities. In this formula- 
tion charges would acquire a Dirac string of electric flux, and monopoles would be 
pointlike. 

Due to Dirac quantization condition, the weak coupling regime with respect to 
charge corresponds to large values of M, or to strong coupling in the dual language 
and viceversa. This is typical of systems with topological excitations. 

The prototype model is the 2d Ising model. The field variable is cr(i) = ±1 and 
the action 

/0jX>(iMi + £) 

Looking at it as a 1 + 1 dimensional field theory configurations like the one in fig. 1 are 
kinks and have a non trivial topology. One can define a dual lattice, by associating 
a point i' of it to each link of the original lattice and a field a*(i') = ±1, assuming 
the value +1 if the sites at the ends of the link have opposite sign, and —1 if they 
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have the same sign. For a kink like the one in figure a* is +1 at the position of the 
kink, —1 everywhere else. 



t 



.............. 

j y { y { y y 



X 

Fig.l A kink in 2d Ising model. 

Duality in this case 15 ) means that the partition function of the dual system has the 
same form of the original one, 

K[a,P\=K[a*,P*] (3-9) 

with the change (5 — »■ (5* ~ 1/(3. 

The two systems are identical, but the low temperature (weak coupling) regime 
of one of them corresponds to the high temperature (strong coupling) of its dual. In 
the ordered phase (a) / 0, (a*) = 0, and in the disordered phase (a*) / 0, (a) = 0. 
(a*) is called a disorder parameter. 

3.2. Monopoles in non abelian gauge theories. 

Monopoles can exist as solitons, or static solutions with finite energy, in non 
abelian gauge theories. They where discovered 16 )' 17 ) in the Georgi-Glashow model, 
a gauge theory with gauge group SO (3) coupled to a triplet of scalar fields (ft: the 
simplest generalization of the U (1) Higgs model. The lagrangean is 

C = —G^G^ + (D^(D^) - V{$) (3-10) 
with the covariant derivative 

Dj> = (d^-gA^ 

and V(4>) the potential 

In the spontaneously broken phase, > 0, (fro = |0o| = /t/. The ansatz 

-> -» r 

do4> = 4>{r) = f(r)<p f f = - (3T1) 

r 

A o (r) = A*(r) = h(r)e iab ^ (3-12) 

gr l 



8 



A. Di Giacomo 



brings to a solution with finite energy, with 

h{r) - =1 f(r) - =1 



The configuration is called an hedgehog, due to the form (3T1) of the Higgs field. 
That this configuration is a monople can be seen by the following arguments. 

a) At large distances, where h(r) = f(r) = 1 a gauge transformation which brings 
cj) on some direction, say the 3 axis, (unitary gauge) transforms the gauge field 
to an abelian field, parallel to 4>, with a Dirac string in space along the 3 axis. 

b) A gauge invariant "electromagnetic field" can be defined 



T v » = 4>- G^u - -4> (D^ A Dj) (3-13) 



9 

Here cf) = 4>/\4>\. 

Ffiu, computed on the monopole configuration, is 

-» 1 f 

(E)i = Toi = H = — k + Dirac string 

gr 2 

H is a colour singlet, and such are the magnetic charges. 

The gauge transformation bringing to the unitary gauge is called abelian pro- 
jection. For the monopole configuration it has a singularity at f = 0, where cf> is not 
defined. 

An alternative way to look at the problem is to use the Body Fixed Frame 
(BFF) lg ) . Usually the same reference frame for colour is used in all points of space 
time, £qj with £q£o = ^ ■> £o ^ £o = £o- Instead three orthonormal unit vectors £i(x) 
can be defined ( with = <5 y , C A£> = I*) with £^(x) = <j>(x). The choice of 
£,2(x) is arbitrary by an angle. The two frames are related by a transformation of 
SO(3), R(x) 

gi{x)=R{x)g (3-14) 
Since (£j) 2 = 1, d^i is orthogonal to £j and 

d„£i = ^A£ i (3-15) 

or 

{&„ - = D^i = (3-16) 

Indeed the body fixed frame changes with x by a parallel transport. 
Eq.(3T6) also implies 

[D li ,D v ]i i = (3-17) 
From the completeness of [D^, D u ] = 0. This means 

G^ v {u) = d jJi u u -d v u IJi + u IJi /\uj v = $ (3-18) 
ui^ is a pure gauge, at least in the regions of space where R(x) is not singular. 



Monopole condensation and colour confinement 



9 



The solution of Eq.(3T6) is 



S(x) = £ 3 (s) = P exp 



Cu^-Tdx^ 1 



£o (3-19) 



which is independent of the choice of the path C if G^ u = 0. P means path ordering 
and T are the generators of SO (3) group, {T 3 )ik = —itijk- 

Expressing £i(x) in terms of the polar angles 9, tp with respect to with polar 
axis along one easily finds 

/ sm9(x)d^ip(x) \ 
^ = -d^9{x) (3-20) 

\- cos9(x)d^(x) J 

At 9 = 0, it, ip is not defined and a singular part of u)^ develops 

to 3 ™ 9 = ( \ (3-21) 
\±d fM ip sin 9(x)J 

and with it a singular field strength tensor 

F^H = ±(d^d u - d v d^ sing Ux) (3-22) 

F^ u is abelian and parallel to (f) = £3. 

For the static monopoles is a Dirac string along the 3 axis at all times: 

the end point of the string is the location of the monopole, i.e. the zero of <j){x). 

The monopole is not an artefact of the abelian projection, but a topological 
feature of the 4> field configuration. 

The monopole configuration described above is also a soliton, and behaves like 
a particle. 

However singularities in can exist also in the unbroken phase of the theory, 
at the zeros of <p, and they are monopoles. Again they are topological properties of 
the <fi field configurations. 

Under infinitesimal gauge transformations exp(iA(x)T) 

w (1 ^w (J + AAw /1 + <9 m A 
A„ -► A„ + A A \ - id M A 

so that 

u» + = gZ» (3-23) 

is covariant. 

In the abelian projected gauge lj^ = and = Z^. The field strength tensor 
can be computed, obtaining 

G^{Z) = G^ U {A) + ~ g FA") (3-24) 
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The gauge transformation which operates the abelian projection is singular, and that 
produces the additional term in G^ v . Also 

Dn(A)$ = (u^ + gAp) A<p = gZ^Aip 

and as a consequence 

-0(D^ A D„<p) = g(Z^ A Z v ) ■ <p (3-26) 
Thus the t'Hooft e.m. gauge invariant tensor 

F^ = $- Gp, - -4> (d^ A D v $) 
reads in the abelian projected gauge 

= d„Zl - d v Zl (3-27) 

since the non abelian term in <\> • G^ u is canceled by the additional term by virtue of 
eq.(3-26). 

The gauge invariant tensor which describes the U(l) field coupled to monopole 
charge coincides with the abelian field of the residual U(l) after abelian projection. 
Because of that sometimes instead of associating monopoles to the field </>, people 
associate them to the gauge which puts <p along the 3 axis. 

Since in QCD the possible fields ip which could define monopoles are many, and 
to each of them a different gauge transformation is associated as abelian projection, 
the statement could be made that monopoles are gauge dependent objects or gauge 
artefacts. This is a misuse of language. 

The physical problem is to identify what 4>(x) fields are relevant to confinement, 
in that monopoles associated to them condense and produce dual superconductivity. 
Monopole species depend on the choice of <$>. However the monopole charges and 
electromagnetic field associated to each (j> are gauge invariant concepts. 

§4. A disorder parameter for dual superconductivity: compact £7(1) 

gauge theory. 

In any field theory in which non trivial topological objects T exist, like monopoles 
or vortices, with a conserved topological charge, a creation operator can be defined 
for them. 

The original construction goes back to ref. 15 ) , and has been developed in different 
forms by many authors 20 ) " 24 \ The basic idea is translation, in the sense of the 
elementary formula 

e ipa \x) = \x + a) (4-1) 

If $(x) is the field describing the theory and < Pt(x,u) is the field configuration pro- 
duced by,say, a monopole sitting at y, then the operator 

fi(y, t) = exp ( i [ d 3 x II$(x, t)0 T (x, y)) (4-2) 



(3-25) 
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with I7$(x,t) the conjugate momentum to 0(x,t), creates a monopole at the site y 
and time t. 

Indeed in the Schrodinger representation \&>(x)) 

li(y,t)\$(x)) = \${x) + $ T {x,y)) 

We will focus on compact U(l) gauge theory in 4d on lattice 25 **' 26 \ where monopoles 
exist and magnetic charge is conserved. 

The building block of the theory is the parallel transport along links of the 
lattice, exiting from site n in direction /j, 

Un(ri) = exp(zea^(n)) = exp(i^(n)) (4-3) 

The parallel transport along the plaquette, the elementary square in the plane \iv, 
is then 

TI^ = exp(i9^ u (n)) 
6^(n) = ApB v - A v 6^ ~ a 2 ef^ (4-4) 
The generating functional of the theory or partition function is 

m-fU^)***-*) (4-5) 

We will choose for S the Wilson action 

S = f3 J2 (1-cosVW) ( 4 - 6 ) 

n,fi<nu 

As P — ■> oo small values of 6^ v are important and 

S^E^ = ^^ 2 C (4-7) 

which is the action for free photons if the identification is made (3 = 1/e 2 . Com- 
pactness of the theory, i.e. the fact that angles 9^{n) only appear as arguments of 
periodic functions in the action, makes Z(f3) and correlations functions of compact 
field variables invariant under the change 

with arbitrary /^(n). A special case are gauge transformations / M (n) = A^<P. 
A critical (3 C exists in the model, (3 C ~ 1.0116 27 )> 28 ). 

For p > [3 C the theory describes free photons. For P < p c electric charge is 
confined. Wilson loops obey the area law, and dual Meissner effect is observed. 
Monopoles exist in this theory. Indeed, since 

from eq.(4-4) 

-4tt < 6^{n) < 4tt (4-8) 
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Since in the plaquette integer multiples of 2ir are not visible, because of compactness, 

^ = 2vrn^ + 9^ - tt < 9^ < tt (4-9) 



one can redifine 9^ v as 



6^ u is the visible flux. 
Now 

obeys the Bianchi identities A^*^ = as can be trivially checked. The visible field, 
however, is 9^ u , because of compactness and can violate Bianchi identities because 
units of 2-7T can be formed which becoomes invisible. A monopole current can be 
defined as 

pff = -^etivpaAunpa (4-10) 

and 



Monopoles are identified as Dirac strings. Since pff is identically conserved, (see 
eq.(4-10)) strings are closed. 

In the pioneering work of ref.( 25 )) such monopoles where numerically detected. 
Their density being higher in the confined phase, (3 < f3 c and dropping to zero 
above /3 C , the density of monopoles was called an order parameter for the transition. 
Although the observed correlation between density and phase is phenomenological 
significant, a genuine disorder parameter for the system should be related to the 
symmetry of the ground state and should then be the vev of a magnetically charged 
operator. A candidate disorder parameter is the vacuum expectation value of the 
creation operator of a monopole. 

In the continuum the general rule eq.(4-2) gives for that operator 



fi(y,t) = exp 



i J (fx E(x, t)-b(x - y) 



(4-12) 



E(x,t) being the conjugate momentum to the field A(x,t), and l/eb the classical 
configuration corresponding to a monopole sitting at y 

Amon{x - y) = -b(x - y) (4-13) 

e 

the factor 1/e coming from Dirac quantization condition is explicitely exposed. With 
some choice of the classical gauge, putting the string along the unit vector n 

l r 27rm r A n 

-b(r) = ^— 4-14 

e e r(r — r ■ n) 

A change of the classical gauge is reabsorbed in the definition (4-12) if E obeys 
Gauss's law. 

On the lattice 

n { = -lmn 0i = -sm6 0i 

e e 
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My, n o) = exp 



p^b^n - y)sm6 i0 (n,n ) 



(4-15) 



The factor j3 comes from the 1/e in the magnetic charge times 1/e in the normaliza- 
tion of the plaquette to the electric field. 

The above definition can be adapted to compactness, i.e. to give a shift of the 
angle instead of its sinus, as follows 19 ^ 



»(y,n ) = exp [s(e l0 (n ) + b\n - y)) - S(e i0 (n )) 
S is the sum of the density of action on the time slice no 



(4-16) 



(4-17) 



2-t e 2 



Since in the limit a — ► 0, C 

a— >u 



(4-18) 



and n coincides with the naive definition modulo a constant coming from the last 
term. 

If more monopoles or antimonopoles are created at time no, b % should be replaced 
by the classical field describing their configuration. 

Correlation functions of monopoles and or monopoles antimonopoles can then 
be constructed. 

We will focus on the correlation 19 ) 



V(x ) = (/Z(0,x o )M0,0)) 



(4-19) 



between a monopole sitting at in space at time 0, and an antimonopole at and 
time xq (propagator of monopole field). 

At large xq we expect, by cluster property 



V(x Q ) ~ Aexp(-M\x \) + (fi} 2 
Translation and C invariance make 

</z(o,z )) = <Mo,o)) = </x) 



(4-20) 



(4-21) 



(/i) is our disorder parameter: (fi) / signals spontaneous breaking of magnetic 
U{1 ) and hence dual superconductivity, as discussed in sect. 2. 

The other important quantity in eq.(4-20) is M, which is the lowest mass of 
excitations carrying monopole charge. The effective scalar field producing dual su- 
perconductivity has a mass larger or equal to M, and hence knowledge of M is an 
important information to determine the type of dual superconductivity. 
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Before going to numerical results we will clarify our definition of fi, by analyzing 
in detail V{xq). According to that definition eq.(4-20) 

V(x ) = ^Z[S + AS] (4-22) 

The factor 1/Z[S] comes from the averaging procedure, Z[S + AS] is nothing but 
the partition function with a modified action: using the definition eq.(4T6) the 
modification consists in replacing 

5(0 <O (O)) S[9 i0 (0) + b\y)] (4-23) 
S(6 l0 (n )) -+ S[e i0 (n ) - b\y)] (4-24) 

Since 

e i0 (fi, 0) = -e\n, 1) + 0*(n, 0) + #°(n + i, 0) - 9°(n, 0) 

the change implied by eq.(4-23) can be reabsorbed by a change of integration variables 

0\n, 1) -»• 0\n, 1) + b\n) (4-25) 

which leaves Z unchanged, because of compactness. 

The result of this change is that 9 l0 (0) is restored to the form it has in Z[S], but 
at no = 1 

9 ij (n, 1) 0^'(n, 1) + Aibj - Ajbi 

A monopole field is added at no = 1, in a form which is independent of the gauge 
choice for 6j. 9 l (n, 1) also appears in 6 t0 (n, 1) and, the change of variables gives 

io (n,l) io (n,l) + bi(n) (4-26) 

which is the same as (4-23), at no = 1. We can now repeat the procedure, and the 
result will be a monopole at time n = 2 and again a change of the form (4-26) 
at time 2. The procedure ends at no — 1 where the change (4-25) is reabsorbed by 
the antimonopole, eq.(4-24). Our construction really produces a monopole at site 0, 
propagating from to x°. 

A direct determination of \i from eq.(4-20) at large value of xq is shown in fig.2. 
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Fig. 2 The disorder parameter vs /3 on a 10 4 lattice. 
Instead of X>(xo) itself it provides numerically convenient to compute 

p(x ) = ^lnP(x ) (4-27) 

or, by eq.(4-22) 

p = (S) s - (S + AS)s+as 

The subscript of the brakets denotes the action used in performing average. Again 
as |xo| — > oo 

p(x ) ~ 2-^ln( / u) + Cexp(-M|x |) (4-28) 

|xo|^oo up 

The typical behaviour of p(xo) near the transition is shown in fig. 3. The typical 
correlation length is of the order of the lattice spacing and M can be determined. 

The behaviour of p^ is shown in fig.4. The sharp drop of (p) around (3 C is 
reflected in a narrow negative peak in p. 
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Fig.3 p(xq) vs xq at (3 = 1.099. The correlation length is of the order of lattice spacing. 
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Fig. 4 p as a function of f3. The sharp negative peak indicates the phase transition. 

For (3 > (3 C the system describes free photons, the Feynman integral is gaussian 
and p can be explicitely computed. Numerically one finds for a lattice L 3 x 2L 



-10.1 L + 9.542 



(4-29) 
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poo tends to — oo as the volume goes large (thermodynamical limit), and 



/ P(a 
Jo 



tends to zero. Above (3 C magnetic U(l) is not broken and (p) = 0. Notice that 
this can only be true in the thermodynamical limit 15 ). In a finite volume (p) is an 
analitic function of (3 and cannot be identically zero on a line of the complex plane 
without being zero everywhere. Only as L — > oo Lee- Yang singularities develop and 
(p) can be zero. 

For (5 < j3 c poo tends to a finite value, compatible with zero as L — > oo, and 
hence (p) / and the system is a dual superconductor. 

The behaviour of (p) around (3 C can be explored by a finite size scaling analysis. 
At P ~ P c a weak first order or second order phase transition takes place. The order 
is controversial 27 -*' 28 ^. In any case the correlation length £ goes large in a range of 
values around /3 C , and an effective critical index v can be defined 



By dimensional analysis 



Z^-APc-P)-" (4-30) 



as ^ — > Zarge the dependence on a/^ can be neglected and 

(p)=p(^0) = f(L^((3 c -(3)) (4-31) 

By use of eq.(4-30) the variable L/^ has been traded with L l l u {f3 c — (3). 
The following scaling law follows for p = ^ ln(/x) 

t = ${l}l"(p c - p)) (4-32) 



This scaling can be matched by appropriate values of i/ and /3 C . We find by best fit 

0c = 1.01160(5) v = 0.29(2) (4-33) 



(3 C agrees with determinations based on completely different methods 27 ) . The quality 
of scaling is shown in fig. 5. 
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Fig. 5 Finite size scaling of p, at the optimal values of is, (3 C . 

We can also compute the critical index by which \x — > at [3 C (p) ~ (/3 C — (3) & finding 

5 = 1.1 ±0.2 (4-34) 

A first order transition would require v = 1/d = 0.25. Further investigation is on the 
way to test if the observed value of v is a finite size effect or it is a true determination 
and the transition is 2nd order. 

We have also measured the penetration depth of the field E, or the mass m of 
the photon. The result is shown in fig. 6, where m is compared to M. 

Further work is necessary to reduce the errors in M. There seems however to 
be sufficient information, at least in the region near /3 C , to conclude that the ratio 
M/m > i.e. that the superconductor is type II. 

Similar results was found by direct observation of the London current in a flux 
tube, in ref.( 29 )). 

If the transition is first order the lattice model does not define a field theory as 
(3 — ► /3 C and the only fixed point is the trivial point f5 = oo, which describes free 
photons. If this is the case changing the lagrangian from Wilson to an alternative 
form, say Villain, also changes the physics of the system, since there is no universality 
class. 

For Villain action the duality transformation can be performed, and condensa- 
tion of monopoles below the critical value has been proved 23 ). In that case it can 
be shown as a theorem that our disorder parameter is equal to that of ref . ( 23 ^ ) , 
although the construction is completely different. Both constructions are based on 
eq.(4-2). The proof of ref.( 23 )) has been extended to Wilson action in ref.( 30 )). 

In conclusion we have a reliable tool to detect dual superconductivity. 
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Fig.6 m (circles) and M (squares) vs (3. 

§5. The XY model in 3d. 

To further check our disorder parameter we have studied by the same technique 
used for U(l) the XY model in 3d, which has a second order phase transition be- 
longing to the same class of universality as the transition to superfluid He^. 

The field variable is an angle 9(i) associated to each site. The action is 

S = /3^^[l-co S (VW)] 

i M 

and the partition function 

z =/n||«p(-5) (5-2) 

The model is compact, so that any change 9{i) — > 9(i) + f(i) with arbitrary f(i) 
leaves correlation functions of compact observables invariant. 

As (3 — > oo S ~ ^(A^6) 2 and the model describes free massless particles. 

At P c ~ 0.454 a second order phase transition takes place, and at (3 < (3 C vortices 
condense. Like in U{\) condensation has been demonstrated in the literature by a 
sharp change of density of vortices 32 ). We shall show instead that for (3 < (3 C the 
U(l) symmetry related to conservation of vortices is spontaneously broken, and we 
will construct a disorder parameter to detect the change of symmetry. 

Let us define 

= d^e (5-3) 



(5-1) 
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The invariance under 6 — ► + / is a gauge invariance since on the transformation 
is 

A lt ^A lt + d^f 

A^ is a gradient, eq.(5-3), and hence a pure gauge 

= 9 n A « ~ = (5-4) 

apart from singularities. 

9{x) can be written as a parallel transport 

6(x) = exp(ieJ A^dx^) (5-5) 

and if eq.(5-4) holds this definition is independent of the path C. 
To investigate possible singularities consider the current 

U = (5-6) 

which is dual to the field strength tensor, is identically conserved (Bianchi iden- 
tity) 

3*% = 

The corresponding charge is 

Q = J d 2 xj (x,t) (5-7) 

or 

Q = J d 2 x (V A A) = <j> Adx = 2im (5-8) 

The path on which the line integral of A is computed is a circle at infinity: the value 
2irn comes from the definition of A„, eq.(5-3). 

In the absence of singularities = and Q = 0. There exist, however, configu- 
rations with Q 0, which are vortices. An example is 

9 q (x -y)=q atan ^"^ 2 (5-9) 

For these configurations 

A = A = t^-^vq 

\x - y\ 

where vq is the unit vector tangent to the circle \x — y\ = r. If the field A is the field 
of velocities the configuration is a vortex with 

Q = <j> Adx = 2-nq (5-10) 

Topology is non trivial. 

As a disorder parameter we will use the v.e.v. of the operator which creates a 
vortex, following the general rule eq.(4-2). Since the conjugate momentum to as 
given by the action is 77 = (3 sin 8q9 the naive definition of \i would be 



H(y, t) = exp 



-(3j2sm(A e(n,t))e g (n-y) 



(5-11) 
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The compact version is 

M (y,t) = exp [s(A 9(t) - ~9 q ) - S(A 9(t)) 

with S the integral of the lagrangean on time slice t 

S(t)=Y / S(n,t) 



(5-12) 



Again, as for the U(l), we compute the correlator 

V(x ) = (fi(6,xo)n(6,0)) ^ ^exp(-M|2;o|) + ^) 5 

|rro|— >oo 



(5-13) 



(//) is the disorder parameter which signals condensation of vortices, i.e. spontaneous 
breaking of the U(l) symmetry. We have from the definition (5-12) 



V{x Q ) 



Z[S + AS] 
Z[S] 



(5-14) 



where S + AS is obtained from S by the replacements at time slices and xq 
respectively 



A 9(n, 0) A 6(n, 0) - 6 q (n - y) 
A 6(n, 0) -» A e(n, x Q ) + 9 q (n - y) 

The change (5-15) can be reabsorbed by a change of variables 

9(n,l)^9(n,l) + 9 q 



(5-15) 
(5-16) 

(5-17) 



which leaves the measure unchanged. However it changes the space derivatives at 
time 1 

Ai9{n, 1) -► Ai9(n, 1) + AiO q (5-18) 



or 



A(n,l) 



(5-19) 



A(n,l)+Ai9 q 

thus adding a vortex to the configuration at uq = 1. 

The other change in the action coming from the change of variables is 

A 9(n,l)^A 9{n,l)-9 q (n-y) (5-20) 

which is like eq.(5T5) with time replaced by time 1. The change of variables can 
be repeated, producing a vortex at no = 2 and a shift of Ao9(n,2) and so on till 
no = xo — 1 when the shift — 9 q of AS disappears with +9 q of eq.(5T6). 
Again it is numerically convenient to study 



p(x ) = ^-V(xo) p p + Ce~ M ^\ 



with 



(fi) = exp 



1 



p{(3')d(3' 



(5-21) 
(5-22) 
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p as a function of f3 is shown in fig. 7 for different lattice sizes. 
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Fig. 7 p vs (3 for different lattices. The peak signals the phase transition. The lines 
at high (5 are the comparison to perturbation theory, eq.(5-23). 

A huge negative peak signals the phase transition. For (3 > f3 c the system describes 
free particles, the Feynman integral is gaussian and p can be computed giving 



-11.33L + 72.7 



(5-23) 



asL^oo,/?^— oo, or (p) — ► 0, as expected in the ordered phase in the thermody- 
namical limit. The agreement with the prediction (5-23) at large j3 is shown in the 
figure: the value (5-23) is represented by the dotted lines. 

For (5 < (3 C p tends, as L — > oo, to a finite value compatible with zero, so that 
(p) / 0: vortices do condense. A finite size scaling analysis around f3 c gives the 
scaling law 



f(L^(/3 c -f3)) 



(5-24) 



v and (3 C can be adjusted to satisfy it. The quality of scaling is shown in fig.8 
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Fig. 8 Finite size scaling for p in 3d XY model. 



We get 

v = 0.669 ±0.065[0.670(7)] (5-25) 
Pc = 0.4538 ± 0.0003 [0.4542(2)] (5-26) 

The values in square parantheses have been determined by completely different meth- 
ods in the literature 33 ) , and are the critical indices of the transition to superfluid. 
The index 5 of the order parameter 

M^iPc-P)* (5-27) 

is also determined 

5 = 0.740 ± 0.029 (5-28) 

Remarkable is the similarity of the shape of p in fig. 7 and fig.4 in spite of the fact 
that the systems have nothing to do with each other. 
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Again we have checked that the construction works, and we will use it to study 
the deconfining phase transition as a possible transition from dual superconductor 
to normal. 

§6. Monopole condensation vs confinement in QCD. 

Contrary to the Georgi-Glashow model QCD has no fundamental Higgs field. 
However, as shown in sect. 3 any local operator <P(x) in the adjoint representation 
has monopoles associated with it, which are located at the zeros of <P(x) for SU(2). 
For gauge groups of higher rank, say SU(3), @(x) is written as a matrix in the 
fundamental representation 

<P(x)=J2® a ^ a (6-1) 

a 

with A a the generators in that representation, monopoles will be located at the sites 
where two eigenvalues of the matrix $>(x) coincide, and will be identified by integer 
charges, in one to one correspondence with diagonal matrices of the algebra with 
integer or zero matrix elements 8 )' 9 \ For SU(N) this means N — 1 U(l) conserved 
monopole charges which can condense in the vacuum. For SU (2) there is one charge. 
The treatment of higher groups does not add any conceptual point, but only formal 
complications. We will therefore use SU (2) formulae to present our arguments. 

The role of <P(x), the operator which identifies monopoles, can be played a priori 
by infinitely many composite operators of the theory: actually by a functional infinity 
of them. Each of them defines monopoles, which can in principle condense and 
produce dual superconductivity. It is not understood a priori to our best knowledge, 
if all of these monopole species are really independent of each other and if many 
of them could condense at the same time in connection with confinement. It could 
also be that U(l) superconductivity for many monopole species is a manifestation 
of a more clever mechanism, a really non abelian superconductivity. This would 
implement the guess by t'Hooft, that all monopole species, defined by any operator 
<P(x) are physically equivalent 10 ) . 

What we can presently do is to investigate these issues on the lattice. Looking 
at the problem from the point of view of symmetry is the most direct way. 

The choices which have been suggested in the literature for <&{x) are 10 - ) 

1) The Polyakov line, i.e. the parallel transport along the time axis to +oo and 
back from — oo via periodic boundary conditions. 

2) Any component of the field strength. 

3) F^Ffj,,,, at least for SU(2>) when it has an octet part. For SU(2) it is a singlet. 

4) The operator which is implicitely defined by the maximization of the quantity, 
with respect to gauge variations 10 )' 34 ) 

Tr {<7 3 fit/ M (n)f2 t <7 3 fi t C/J(n)fi} = max (6-2) 

with Q(n) a generic gauge transformation. This procedure defines an operator 
<P in the adjoint representation which coincides with 

J2 (u^(n)azU^n) + U^n - A)cr 3 C/J(n - n)) (6-3) 
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in the maximal abelian gauge. The explicit form in a generic representation is 
not known. 

To explore dual superconductivity we will construct a disorder parameter, as the 
v.e.v. (fi) of an operator fx with U(l) magnetic charge, as we have done for U(l) 
and XY model. 

Whatever the choice of <P, in the gauge in which <P ■ a is diagonal, i.e. after 
abelian projection, any link U^{n) can be written 34 ) 

e io- 3 o M (n) e «o- 2 7 M (n) g ~ i(J3a M (n) 

e i CT3 e M (n) = e ja 3 (Q M (n)+/3 M (n)) ig the abelian link) which parallel transports the U(l) 
field related to the monopole charges defined by <P. 

The creation operator of a monopole at time t = will be defined by changing 
the kinetic term of the action at that time by adding the field of the monopole to 

0*)( n )35)_ 

The change on the i — plaquette ZZjo will be 

II l0 (n,0) ^ n' t0 (n,o) b 



= io-3(a M (n)+/3 M (n)) 



(6-4) 



il i0 (n, 0) = Tr [Ui(n, 0)U (n + i),0)u}(n, l)[/ f (n, 0) 
n' i0 (n, 0) 6 = Tr f[//(n, 0)[/ (n + »), 0)f//(n, l)[/ f (n, 0) 



and 



C//(n,0) = e iA{n)a3 Ui(n,0)e ib ^ a3 e 



■iA(n)crz 



(6-5) 



We have operated the separation of the vector potential describing the monopole 
into transverse and longitudinal part 

bi{n) = bi{n) + diA(n) dib^n) = 

The gauge part e lA ^ a3 and e ' lA ( n + 1 )' T 3 can b e reabsorbed by a rotation of the C/o's 
which leaves the functional measure invariant, so that the definition is independent 
on the choice of the classical gauge for bi(n), and the net effect is to add bf{n) to 
the abelian phase of Ui(n, 0). 
Consider now the correlator 



As usual we can write 



V{t) = {fi(x,t)ti(x,0)) 
Z[S + AS] 



V{t) 



Z[S] 



(6-6) 
(6-7) 



where S + AS is obtained by the substitution 

il 0i (n,0) 7I fe 0i (n,0) 
n 0i (n,x )^n°\(n,x ) 
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The subscript b, —b is to recall the sign of the monopole charge. 

The effect of this procedure is to have a monopole created at t = 0, which 
propagates to time t, when it disappears. The construction to show this is identical 
to U(l) if <P is the Polyakov line. There after abelian projection the temporal links 
are diagonal 

U (n) = e ia3a °W 

and the operator e %azh ^ in the definition of U- commutes with them. 

Then a change of variables from Ui to U[ which leaves the measure invariant 
brings n i0 back to the original form, but changes iTj^n, 1) by adding a monopole to 
the abelian field 

AiOj - AjQi -► AiOj - AjO, + Aibj - A,b, 

The other change is to ZP°(n, 1) — > i7£°(n, 1). The procedure can then be repeated 
till at time t the change is reabsorbed by —b. 

If the Uq's are not diagonal in the abelian projected gauge, the construction is 
the same modulo additional parallel transports in the change of variables, which do 
not modify its content. 

Again one can either measure V itself or 

p(t) = A in V(t) = (S)s-(S + AS) S+ AS (6-8) 

One expects 

V(t) = (p,) 2 + Ae~ Mt (6-9) 

pit) ~ p + Ce- Mt (6-10) 

P=2^1n( M ) (6-11) 

A similar analysis to the one performed for U(l) brings to the determination of (p) 
in the thermodinamical limit, and of M. (p) ^ signals dual superconductivity. 
A measurement of the penetration depth of the field allows to establish the type of 
superconductor. 

At finite temperature, i.e. keeping the time extension Nt of the lattice much 
smaller than the space extension A^^, the correlation of (p) to confinement can be 
studied. There a single monopole is used, and (p) is measured. A careful examination 
of the construction given above shows that, when computing Z[S + AS] of eq.(6-6) if 
periodic boundary conditions are used, the change of variable described above adds 
more and more monopoles when we go trough the boundary in time. The way to 
have one monopole is to use antiperiodic boundary conditions in time. This suggests 
that monopoles behave as fermions. 

We are systematically exploring by our disorder parameter the deconfining tran- 
sition in SU(2) and SU(3), by different choices for ^, the field which defines monopoles. 
We are measuring penetration depths and critical indices. As for <P we consider 
Polyakov line, field strength component and max abelian projection. Typical be- 
haviour for p are shown in fig.9 and fig. 10 for SU(2) and SU(3) 
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Fig.9 pvs [3 for SU(S) gauge theory. The peak signals deconfining phase transition. 
Here monopoles are defined by the abelian projection on Polyakov line. 

Preliminary evidence is that for all the species of monopoles considered, vacuum 
behaves as a dual superconductor, and undergoes a phase transition to normal at the 
deconfinement point. This supports the guess of t'Hooft about physical equivalence 
of different monopole species 10 ). 
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Fig. 10 Same as Fig.9 for SU{2). 

§7. Concluding remarks. 

Our strategy to answer the question if dual superconductivity is the mechanism 
of colour confinement, is to look at the symmetry of the vacuum. For that we have 
constructed a disorder parameter, which directly detects dual superconductivity 

The construction has been tested in known systems, like the U(l) compact gauge 
theory and the XY model in 3 dimensions. 

We have evidence that for many choices of the effective field <P(x) defining 
monopole species, dual superconductivity is present in the confined phase, and dis- 
appears in the quark gluon phase. 

Additional relevant information from lattice is that 

1) Flux tubes exist in the space between propagating QQ pair 36 )' 37 ). 

2) If one single species of monopoles were at work to produce superconductivity, 
then the electric field in the Abrikosov tubes should be the field of the U(l) 
group to which monopole charges belong. An analysis of the colour content of 
the flux tubes shows instead that its direction in colour space is uncorrelated 
to the direction of S 37 ^ 38 \ 

Moreover, whatever the abelian projection is, there exist one gluon in SU(2), two of 
them in SU(3), which have zero electric charge with respect to the residual £/(l)'s, 
and therefore cannot be confined. The adjoint string tension is zero, and also this 
fact seems to contradict lattice observations 37 ) . 

A different approach to the problem is to pay less attention to symmetry, and 
look at more quantitative facts, like abelian dominance u ) and monopole dominance 
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12) _ 

The abelian part of the field as defined by abelian projection, in the max abelian 
gauge, is a good approximation to full dynamics, and of it the contribution of 
monopoles is dominant. In a sense abelian dominance is expected, since, after max- 
imal abelian projection the links are diagonal within 85%. However the fact that it 
happens is surely relevant. 

Looking at symmetry, as we do, is a complementary approach. 

Hopefully a picture will emerge from all these efforts, which will improve our 
understanding of the theory. 
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